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Abstract. Carbon nanofibers (CNF) are used as components of planar photonic 
crystals. Square and rectangular lattices and random patterns of vertically aligned 
CNF were fabricated and their properties studied using ellipsometry. We show that 
detailed information such as symmetry directions and the band structure of these novel 
materials can be extracted from considerations of the polarization state in the specular 
beam. The refractive index of the individual nanofibers was found to be ucnf = 4.1. 
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Two-dimensional (2D) photonic crystal (PC) slabs are materials with a periodic 
modulation of the refractive index in a 2D plane and a finite extension in the 
perpendicular direction. These systems have generated much interest in recent years 
as a means to control light at length scales on the order of the wavelength [1] . 

In this article we report the study of properties of 2D PC slabs consisting of free- 
standing carbon nanofibers (CNF) grown on a metallic substrate. The optical properties 
have been investigated using the technique of variable-angle spectroscopic ellipsometry 

PCs based on CNFs are thought to be promising systems for the development of 
tunable PCs in the visible range |3], making use of the favorable actuation properties 
of the CNF. In order to fully exploit this potential application, it is important to first 
understand the static optical properties of such 2D arrays. It was earlier demonstrated 
that vertical CNF can act as optical antennas to transmit light in a fashion similar to 
the transmission of radio waves in RF applications Two-dimensional lattices of CNF 
have also been shown previously to act as diffraction gratings P,[T]. However, there has 
been no experimental determination of the band structure of the CNF-based PC and 
corresponding detailed comparison with theoretical predictions. 

The results reported here demonstrate that it is possible to map out large regions 
of the band structure of the CNF-PC using ellipsometry. The measurements are in 
good agreement with finite-difference time domain [5] band structure calculations using 
a refractive index of ticnf = 4.1 for the CNF. It is also shown that the symmetry 
directions of the samples can be easily detected. 

Two-dimensional PC slabs consisting of free-standing CNFs are directly grown on 
a metallic underlayer using plasma-enhanced chemical vapour deposition (PE-CVD). 
The technique, with some modifications to the growth parameters, has been described 
previously [9]. The lattice structure is defined by using electron-beam lithography to 
pattern Ni catalyst dots with the required lattice parameters. Arrays were patterned 
with lattice constants, a, ranging between 500 nm and 200 nm in different configurations 
along the x- and y-directions. Random samples were patterned so that the average CNF 
density corresponded to that of the square lattices. As an example, for comparison with 
the a = 500 nm lattice constant, CNFs are placed randomly with 0.8 % filling factor. 

Titanium substrates coated with a 20 nm titanium nitride layer were used as the 
underlying support. The size of the Ni catalyst dots was 50 x 50 x 30 nm^ {in x x y x z 
directions respectively). The growth chamber was first pumped to 10^^ Torr before 
ramping the temperature by 100 °C/minute while ammonia was introduced to the 
chamber at 60 seem. The DC plasma was ignited at a pressure of 4 Torr and a substrate 
temperature of 500 °C. This step is used to clean the substrate and activate the catalyst 
particles and is carried out for 2 minutes with a plasma current of 4.3 mA cm~^. During 
the activation step the temperature was raised to 600 °C at a rate of 50 °C/min. At 
600 °C, acetylene was introduced to the chamber at 15 seem and the plasma current 
reduced by 33%. The temperature was ramped to 700 °C with a ramping speed of 
50 °C/min. Growth time was adjusted depending on the desired height of the CNFs. A 
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height of approximately 1 /iin is reached after 20 minutes. The CNFs were found to grow 
via a tip-growth mechanism with the Ni particle situated at the tip after growth. The 
samples used for the ellipsometry measurements were ca. 1 fim long, 50 nm diameter 
CNFs. An example of a square lattice and of a random pattern are shown in panels 
(c) and (d) of Figure [U along with a schematic illustration of the PC (b) and a sketch 
illustrating the experimental geometry (a). 

In reflection ellipsometry, the polarization change of light reflected off a sample 
surface is measured. This change depends on the optical properties and the structure 
of the sample, including any overlayers. The basic quantity measured is 

P=t^ = tanvl/e^^ (1) 

where Rp and Rg are the complex- valued reflection coefficients for light polarized parallel 
(p) and perpendicular (s) to the plane of incidence. \E' and A are called the ellipsometric 
angles. The ratio of amplitude change for p- and s-polarization is given by tan\E'. The 
corresponding difference in phase changes is given by A. 




(c) (d) 

Figure 1. (a) Ellipsometry setup with angle of incidence, 9, and sample rotation 
angle, marked, (b) Sample illustration with layer structure. SEM micro graphs 
of (c) square lattice with 500 nm lattice constant (~ 0.8 % filling fraction) and (d) 
random pattern with the same CNF filling factor as the sample in (c). Scale bars in 
(c) and (d) are 1 /j,m. 
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However, if a sample is anisotropic, it is generally not possible to fully characterize 
it with ordinary ellipsometry, which only takes into account the diagonal elements of 
the reflection Jones matrix [10], as was done in an analysis of 2D Si PC by Hsieh et 
al. [11]. When the off-diagonal elements of the Jones matrix are non-zero, generalized 
ellipsometry is needed. This will ensure that conversion between s- and p-polarization 
is correctly taken into account. In generalized ellipsometry, three complex- valued 
parameters are defined as 

tanvl/^,e^^- 

tanvl/^,e^^- (2) 
tanvl/,,e^^-, 

where Rpp and Rss are the reflection coefficients for the p- and s-polarizations, 
respectively, and Rps and Rsp represent the coupling from p- to s-polarization and from 
s- to p-polarization, respectively. In this study only \l/pp data (denoted simply as \& 
for convenience) are presented. This choice is based on the fact that \E' is the most 
fundamental parameter in generalized ellipsometry. It is also the easiest parameter 
to start analyzing and has the closest connection to PC structure. A more detailed 
ellipsometric analysis, considering a larger part of the Mueller matrix, will be presented 
elsewhere. From (Mi we find that 



tan\& 



Rpp 



Rg 

which thus represents the ratio between Rpp and Rss- 

The ellipsometer used in this study is a dual-rotating-compensator ellipsometer 
from J. A. Woollam Co. Inc. The instrument provides all 16 elements of the reflection 
Mueller matrix from which the generalized ellipsometry parameters can be derived [2]. 
The full Mueller matrix was measured at different angles of incidence, 6, and different 
rotations of the sample, u (see Figure [Tj^a) for angle definitions) . 

For the 500 nm x 500 nm lattice, the ellipsometry parameter measured at two 
sample orientations, u = 0° and u = 45° (where i/ = correspond to the x-direction 
of the sample), is plotted against wavelength in Figure |2l Measurements from the bare 
substrate and the random pattern are also shown for comparison. The data from the 
bare substrate show no pronounced features (except for the dip around 300 nm which is 
due to the low refractive index of titanium at that wavelength [13j). The random pattern 
shows a fairly smooth variation in \1/ but with a value that is significantly lower than 
for the bare substrate. This indicates that p-polarization is reflected much less strongly 
than s-polarization for the random pattern compared to the bare substrate. This effect 
can be related to the coupling of the incident p-polarized light to the surface plasmon 
polaritons (SPPs) [8] supported by the metallic substrate. Note that this coupling is 
only possible due to the scattering of the incoming fields with the disordered CNFs at 
the substrate surface. 
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Figure 2. Ellipsometric parameter ^P, see ([T]), plotted against wavelength for the 
square lattice with lattice constant 500 nm, i/ ~ 0° (solid black) and i' — 45° (dashed- 
dotted red), and for the random pattern (dotted green) with the same CNF filling 
factor. Also ^ for the bare substrate is shown (dashed blue) as a reference. All curves 
are shown for angle of incidence 9 — 75°. The second (first) peak that appears for the 
square lattice for = 0° {i^ = 45°) is marked with an arrow and denoted I (II). 



Peaks appear in the \1/ vs. A plot for the square and rectangular lattice structures as 
a consequence of their periodicity. Due to the regular lattice, the system support bound 
electromagnetic modes characterized by their optical band structure [13]. Since the 
CNFs are narrow (compared to the lattice constant) they represent a weak modulation 
of the metallic substrate surface. This implies that the first (lowest) band will not be 
much affected by the CNF modulation, remaining very similar to the dispersion relation 
of the SPPs supported by the bare substrate. However, at the Brillioun zone (BZ) edge 
a photonic band gap will open due to the presence of the CNF PC. 

The group velocity in the PC is related to the modal dispersion via the relation 

Vg(k) = VkW, 

where k denotes the photonic Bloch wavevector in the PC, and u the angular frequency. 
The resulting effective refractive index Ues is direction dependent (a tensor) due to the 
anisotropic dispersion relation. In particular, at the BZ edge dco/dk — ?■ 0, corresponding 
to a diverging Ugs in the substrate plane, which yields an impedance mis-match between 
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the vacuum and the PC regions and causes reflection, mostly for the p-polarized hght 
whose dispersion relation is more strongly affected by the PC. Therefore, we can expect 
sharp spectral features at wavelengths for which the dispersion relation of the modes 
supported by the structure yields fc-vectors at the 2D BZ edge. If we consider the 
data shown in Figure [2] for the square 500 x 500 nm^ lattice, the longest wavelength 
(smallest k) maxima for i/ = 0° and u = 45° appear at ca. 950 nm and ca. 690 
nm, respectively. The corresponding wave vectors, 6.6 /im~^ and 9.2 fim~^, are simply 
related by cos (45°) ^ 0.71 as one would expect from considering the geometry of the 
BZ. 

It should be noted that it is the in-plane refractive index that diverges. This can, 
however, still give rise to an increase in reflection in the 2;-direction, see e.g. |15j . 

A rectangular 2D PC has a reciprocal lattice that is also rectangular [16], see 
inset in Figure El Here we consider rectangular lattices with CNF at positions 
Hpq = p 500 nm x + q 400 nm y where p and q are integers. From such geometry, 
we expect the flrst spectral maximum to correspond to an edge of the BZ that is 
either parallel with the y-axis (perpendicular to direction (10)) or with the x-axis 
(perpendicular to direction (01)). When the azimuthal angle u is varied, we expect 
the peak position to shift according to 

,(0) 



''peak 



^(10) - 500 nm 



tanz/ < 



cosu 400 nm 

^(01) 500 nm 

^ — -, tanz/ > . 

sin u 400 nm 



Here k^^^-^ is the wave vector of the first peak for a rotation angle u such that the specular 
beam in the ellipsometer is in the (lO)-direction. The parameter k^^\^ is similarly 
obtained for a specular beam in the (01) direction. From the geometry of the BZ it 
follows that 

L.m _ 5 (0) 

(01) ~ 4 (10)- y"^) 

In Figure [3l the lowest wavelength peak is plotted together with the predicted fcpeak- 
Remarkably, the agreement between the peak position and this prediction is excellent. 

Next, we wish to understand the physical origin of the other peaks that appear in the 
spectra. The band structure is calculated using the finite difference time domain method 
[T7] . Although the geometric structure of the PC is well defined (Figured]), the refractive 
index of the CNFs is not known. We therefore calculate the band structure for several 
different values of the refractive index and compare the results to the ellipsometric 
data. Best agreement is found for a choice of ncNF = 4.1 0. The band structure for 
a 2D square lattice PC calculated using this value of ncNF is shown in Figure HJ Note 
that for PC with high-index pillars embedded in low-index material s-polarization in 

I This can only be considered a phenomenological model. For instance, the refractive index of the 
CNF cannot be expected to be constant over such a wide frequency range and no dissipative effects 
associated with the conductivity of the CNF have been included. 
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Figure 3. Experimental (circles) and calculated (curve) fc-vector magnitude of the 
peak maximum as a function of the azimuthal angle, for the first peak (longest 
wavelength) in the 500 x 400 nm^ lattice, corresponding to the longest wavelength 
peak seen in Figure [2] for the square lattice. 

the ellipsometry setup is only weakly affected while the p-polarization in ellipsometry is 
strongly affected [T8] . 

As can be seen in Figure HJ for p-polarization, there are several regions of {u, k) 
space where the band structure is flat, leading to a divergence in the effective refractive 
index, ries- The high effective refractive index causes strong reflection of the p-polarized 
light, increasing Rp in ([1]), and hence showing up as peaks in \E'(A). Since the s- 
polarization is only weakly affected by the CNFs, we basically see only peaks and no 
dips in as were observed for the Si PC reported by Hsieh et al. [H]. 

To test if the detected peaks in \E' can be related to the band structure, was 
recorded for an angle of incidence 6 = 75°. The in-plane wave vector k\\ was obtained as 
k\\ = sin 6 k where k = u/c is the free space wave vector. In Figure HI the experimental 
data are shown as empty red dots superimposed on the calculated band structure (solid 
green dots). The figure demonstrates how our experimental technique allows us to map 
the dispersion relation of the electromagnetic modes supported by the CNF PC. Note 
that the aim of this work is to test the capability of ellipsometry measurements to study 
guided modes in 2D PCs and that a more exhaustive reproduction of the theoretical 
bands is out of the scope of this paper. Our experimental results are in good agreement 
with theoretical data even in the XM axis, which is harder to probe experimentally. 

The theoretical bands within the shaded region in Figure H] correspond to guided 
modes that lie below the lowest light line. These modes do not experience any radiation 
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Figure 4. Band structure of a 2D square PC consisting of dielectric pillars with radius 
of 0.05a and a refractive index of 4.1, calculated using finite difference time domain 
computations. The circles mark peak max of 4* for different rotations. The blue dashed 
circles I and II correspond to peaks I and II in Figure [2j Inset shows the BZ of the 
square lattice with high symmetry points marked and the irreducible BZ highlighted. 
Band structure is filled green circles and solid black lines are light lines. 



leakage while propagating along the PC structure, which implies that, in principle, they 
cannot be excited by ellipsometry measurements. Note that energy and momentum 
must be conserved in the scattering of the incident light with the sample. This explains 
why these modes remain hidden in the measurements along the F-X and M-F. Along the 
X-M direction, however, our experiments seem to reproduce the second lowest dispersion 
band (below the light line). The appearance of non- leaky modes supported by 2D PCs 
in reflection experiments has been recently reported by Paraire and co-workers [I2]. In 
this work, this effect is linked to the presence of imperfections and to the finite size of 
the experimental sample. 

The solid black lines in Figure H] indicate the light lines. These correspond to 
the dispersion relation of grazing Bragg waves that propagates parallel to the sample 
surface. Note that, although these modes are not bounded to the PC structure, they play 
a relevant role in the scattering of the incoming light with the sample in our ellipsometry 
experiments [18] . It should be noted that the guided PC modes tend to the light lines 
and eventually overlap at certain wavelengths within the BZ (see Figure H]). 
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In Figure [2], two large peaks are marked as I and II. These are depicted in Figure H] 
as ellipses. They show up for the PC rotated in high symmetry directions, corresponding 
to directions with strong dispersion, as can also be seen in Figure HI Hence we conclude 
that ellipsometry also can be used to map symmetry directions of PCs when these are 
unknown. 

In conclusion, there are no complete band gaps for any polarization in our samples 
due to the low filling factors of the PC under study. However, we still see strong effects 
at positions where the bands flatten out. The measured effects are much stronger than 
one might initially expect but they are in agreement with earlier computational results 

Ellipsometry is a non-destructive, well-established technique for characterizing thin 
films. In this paper we have shown that it is possible to use ellipsometry to characterize 
2D PC slabs and determine their band structure. It is also possible to easily identify 
symmetry directions and lattice parameters. 
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